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ABSTRACT 



We use the formalism of lSzapudl lj2004f ) to extend the theory of large angle redshift 
distortions. We identify a previously neglected term in the Jacobian, which contains 
a purely geometric non-linearity while still being linear in terms of the dynamics. We 
perform linear perturbation theory while keeping these non-perturbative terms intact. 
The results represent a significant correction to previous calculations and are in excel- 
lent agreement with our measurements in the Hubble volume simulation. It appears 
that the new formulae amount to a satisfactory description of redshift distortions on 
large scales. 



1 INTRODUCTION 



Redshift distortions represent a curse disguised as a blessing 
for high precision cosmological applications. Radial coordi- 
nates of redshift surveys contain limited phase space infor- 
mation, which in principle can be used to constrain theories 
more than configuration information alone; moreover, veloc- 
ities are sensitive to structure outside of the survey bound- 
aries which potentially translates into a larger "effective vol- 
ume". On the other hand, redshift distortions are plagued 
with non-linearities, both on large and small scales, there- 
fore in worst case they could amount to poorly understood 
contamination of the configuration space data. Our aim is 
to extend the theory of linear redshift distortions such that 
large angle information could be successfully extracted from 
galaxy surveys. 

The work of lDavis fc Peebles! (|l983l ) and lPeeblej (|l980h 
showing that redshift distortions affect the power spec- 
trum spawned a lot of activity. The all-im portant linear , 
plane-parallel limit was first calculated by iKaiserl (| 1987h . 
showing that the effect on the power spectrum corresponds 
to "squashing". The other well known "fingers of God" 
effect dominates small scales, and is irrelevant for our 
present study. The Kaiser formula has been generalized 



but if our aim is to extract as much information as possible 
from a given su rvey, it woul d be d esirable to add them in. 
This motivated ISzalav et al.l (|l998l ) to perform a full linear 
perturbation theory calculation without the distant observer 
limit. The calculation in this results a finite expression only 
in coordinate space: in Fourier space an infinite series will re- 
sult for the redshift distorted analog of the power spectrum. 
These formulae were later further gener alized to include 
high-z effects in various cosmologies by iMatsubara et all 
(2004). These results provided the theoretical foundation of 
several subsequent analys es of wide angle redshift su rveys, 
such as iPope et ail i|2004l ) and lOkumura et al.l (120071 ). De- 
spite its elegance, t he theory did not ag ree well with dark 
matter simulations. IScoccimarrol (|2004| ) pointed out that 
this might be due to non-perturbative effects. 



for r eal space soon after (e.g., iHamiltonl 1 19931 : ICole et al.l 
1 19951 ). These theories have been used to analyze surveys 
such as the Point Sou rce Catalog Redshift (PSCz) Sur- 
vey (|Tadros et aill999l ) . the Two-Degree Field Galaxy Red- 
shift Survey (2dFGRS: IPeacock et al.ll200ll : lHawkins et ail 
120031 : iTegmark et al.ll2002l). an d the Sloan Digital Sky Sur- 
vey (SDSS: IZehavi et alj|200i ). 

The distant observer approximation only holds if pairs 
are separated by a small angle. This means that a large frac- 
tion of pairs needs to be thrown away from modern wide 
angle redshift surveys when they are analyzed in this limit. 
These pairs are typically fewer and noisier than close pairs, 



ISzapudil l|2004h reanalyzed the redshift distortion prob- 
lem from group theoretical point of view showing that tripo- 
lar spherical harmonics provide an excellent basis for expan- 
sion, and result in especially compact formulae. In addition 
it provided specific coordinat e systems, one of w hich recov- 
ers the Legendre expansion of Sz alav et al.l (| 199Sf ) . while the 
other represents the same information in an even simpler 
Fourier mode expansion. We use this formalism to take into 
account another term in the Jacobian neglected by all pre- 
vious work. This term, while linear in terms of the small 
fluctuations, is essentially non-linear from the point of view 
of geometry: it contains a 1/r prefactor, which was also 
the justification for neglecting it. Moreover, this term, if 
expanded in terms of bipolar spherical harmonics (or any 
other way), would contribute infinite terms. In this paper 
we introduce a hybrid approach, where we leave the essen- 
tially non-perturbative terms in the expansion, i.e. our tripo- 
lar expansion coefficients will still contain angular variables. 
As we show later, this procedure results in a finite num- 
ber of terms, and it provides dramatic improvement in the 
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agreement with simulations. In retrospect, the omission of 
this term, while intuitively reasonable, is not justified, as its 
contribution can be critical on the most interesting scales of 
tens of h _1 Mpc's. 

In the next §2 we present the theory of linear redshift 
distortions including results from the newly identified geo- 
m etric term in th e Jacobian. We follow closely the formalism 
of lSzapudl (|2004T l. mainly focusing on the new aspects of this 
calculation. For reference, we print the full result, which has 
about twice as many terms as previously. In §3 we compare 
our results with preliminary measurements in the Hubble 
volume simulations, and present our conclusions. 
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We use x for denoting x\ — X2- On the right hand side one can 
find the Wigner 3j symbols and we define the normalized 
spherical functions as C; m = yj 4tt/21 + lYi m ; these latter 
result in simpler expressions. 

Eq. ((2| has become more complex with the additions, x\ 
and X2 appear in the denominator resulting in the following 
angular dependence 
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2 REDSHIFT DISTORTION OF THE 

TWO-POINT CORRELATION FUNCTION 

We use linear perturbation theory to predict the redshift dis- 
torted two-point correlation function in terms of the under- 
lying power spectrum. Our calculati on is based dir ectly on 
the tripolar expansion formalism of ISzapudU (|2004T ). there- 
fore our focus will be on the additional terms arising from 
the Jacobian. 

The exact mapping between real and redshift space is 
Si = Xi — fvjXiXj, where the "hat" denotes the proper unit 
vector, / = Q-g— and the velocity has units which provides 
that its divergence is equal to the density up to linear order. 
From this, one can calculate the derivative of this matrix: 
dsi/dxk = Sik + Oik(v) where O is linear in v. This results in 
a linear Jacobian J = 1+TrO = 1 — fxiXjdiVj — 2/^-J 2 -. The 
last term in the previous expression is usually omitted due to 
the fact that it scales with 1/x, i.e. it would tend to zero for 
large distances, which loosely correspond to large angles as 
well. Closer examination of this term shows that it is of the 
same order as the previous term, not only in perturbation 
expansion (linear), but also in order of magnitude. Our goal 
is to propagate this new term through the full calculation. 

The linear density contrast and the two-point function 
can be expressed in the usual fashion. 



Expanding these terms into tripolar spherical harmonics 
would yield infinite terms, but simplification arises from the 
fact that they can be fact ored out of th e integrals. All the 
rest can be expanded as in lSzapudU (|2004T ). resulting in finite 
expressions. We introduce 0i to denote the angle between &i 
and x and 02 for the angle between £2 and x. We emphasize 
that the coefficients of this (quasi-)tripolar expansion still 
has an angular dependence in the form of gi and g2- 



£s = ^ B llhl (x,(f) 1 ,(j>2)Sl 1 l 2 l{x 1 ,X2,x). 
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After performing the expansions, only a finite n umber of 
coeffic ients survive. For reference, the ones from ISzapudil 
d2004h are: 
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where Pi and P2 are Legendre polynomials and P(k) is the 
linear power spectrum. The third term in each of the brack- 
ets correspond to the extension of the previous results; these 
would tend to zero in the plane parallel limit. At wide an- 
gles, the separation between the galaxies and the distance 
between a galaxy and the observer are of the same order, 
therefore kx is of order unity. This shows explicitly that the 
order of this term can be as large as the previous, and the 
detailed calculation confirms this. 

Next we express the angular dependence of the correla- 
tion function with tripolar spherical harmonics. 

5( 1 ! 2 i(xi, X2, &) 



and the new terms, the main result of this paper, are 
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where ^{x) = J dk/2n 2 k m j l (xk)P{k) with j being the 
spherical Bessel function. 
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Figure 1. (Left) The measurement of the correlation function without redshift distortion of the Hubble volume simulation (dots) 
compared with linear theory (dashed and solid lines). The errorbars were estimated from 9 3 subvolumes of th Hubble volume. Shifting 
the theory by 0.001 downward, motivated by the integral constraint, provides an excellent fit to the data. (Right) Redshift distorted 
correlation function of the Hubble volume simulation (dots) at constant opening angle (0.71 radian) and while the ratio of the distances 
of the particles in the pair are kept fixed (at 1.57). The errorbars were estimated as before. The lines indicate the linear theories with 
(this paper) and without the geometric terms. The solid line is the corrected theory with a downshift of 0.0024. The integral constraint 
correction is expected to be larger since the average of the two point function is larger. 



For furth er elabo r ation we choose coordinate sys- 
tem a) from ISzapudil (|2004T ) . This corresponds to our 
previous choice of angles with 0i, 02, with which the 
Si 1 i 2 i{x 1 ,x 2 ,x) = Si lh i(n/2,4>i,n/2,(p2,n/2,0) functions 
can be expresse d using cosines and sines only. Using the 
same notation as ISzapudil l|2004 ): 

£ 8 (0i,02, a:) = a ni „ 2 cos(ni0i) cos(n 2 02) 

n l i 7L 2 —0,1 , 2 

+b ni „ 2 sin(ni0i) sin(n 2 02). (9) 

Again, for reference, the previously calculated coefficients 
are 
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and the new expressions of this work correspond to 
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It is worth to emphasize again that the angular depen- 
dence gi and g2 is suppressed for clarity in the above for- 
mulae, but it is obviously carries through according to the 
definition of these functions. If the equivalence of the con- 
figurations (0i, 02) — > (tt — 2 , n — 0i ) is taken into account 
(same pairs can be counted twice), the number of indepen- 
dent new coefficients is five, i.e. the number of terms ap- 
proximately doubled. Next we explore the relevance of these 
calculations, and compare the theoretical predictions with 
measurements in dark matter only TV-body simulations. 



3 DISCUSSION AND SUMMARY 

To understand our results, we expanded our formulae to 
identify leading order corrections to the Kaiser limit. 

The leading order corrections to the distant observer 
approximation are second order. Using the notation i(0i + 
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terms in A0 results in 



A0, and keeping leading order 
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The first line of eq. (|12[) corresponds to the Kaiser formula 
(Aa = 0).The next line contains leading order corrections 
corresponding to previous work only, and the third line col- 
lects leading order corrections from the new term in the 
Jacobian. These are all of the same order, reassuring the 
need of keeping the geometric non-perturbative terms. We 
conjecture that the terms containing the cot (<f>) could be 
responsible for the reported failure of the linear theory f or 
small angles along the line of sight (lOkumura et al.l l2007). 

As a preliminary test of the validity of our calculations, 
we measur ed correlation func tions in the Hubble volume 
simulation (|Evrard et al.l [20021 ), using cosmological param- 
eters as = 0.9, n s = 1, O m = 0.3, Six = 0.7, h = 0.7, 
n h h 2 = 0.0196 and a volume of (3000Mpc/i -1 ) 3 , with and 
without redshift distortions. The volume of the simulation 
was divided into 9 3 subvolumes to obtain the errorbars. 

The left panel of the figure shows the measured and 
the theoretical two-point functions without redshift distor- 
tions. The theory agrees with the measurements only after a 
shift by a consta nt. This could be due to the "integral con- 
straint" problem ([Peebles 1980, e.g.,), possibly compounded 
with slight non-linear effects. 

Next, an observer was placed at the center of each sub- 
volume and the mapping between real and redshift space 
was performed using the velocities recorded in the simu- 
lation. The correlation function was then measured using 
brute force counting of pairs in high resolution bins match- 
ing our choice of coordinate system described earlier. The 
right panel of the figure presents wide angle redshift dis- 
tortion theory both with and without non-perturbative ge- 
ometric corrections. The theory without corrections cannot 
be made to agree with the measurements even using a con- 
stant offset. In contrast, the theory presented in this paper 
provides excellent agreement with the measurements if the 
effects of integral constraint are taken into account. Note 
that the shift corresponding to the latter is expected to be 
larger with redshift distortions included simply because the 
two-point function is enhanced on large scales. 

While these measurements are preliminary in the sense 
that we did not try to span the full parameter space of wide 
angle redshift distortions, the results presented in this figure 
are typical in the sense that other configurations we mea- 
sured showed similar improvement. Scanning the full pa- 
rameter space with our present brute force two-point cor- 
relation function code would be impractical, since we need 
a very large number of pairs in each bin to beat down the 
errorbars enough that the difference between the two theo- 
ries can be reliably measured. Although we developed a fast 
grid based code as well, we found that at these small values 
of the correlation function the pixel window function effects 
become important. These are more complex for the redshift 
distorted correlation function depending on three variables 



than in real space. Such effect should be modeled very ac- 
curately before one could fully span the available parameter 
space. Such detailed modeling together with any other real- 
world issues affecting cosmological parameter estimation is 
left for subsequent research. 

With these caveats we conclude that our theory of wide 
angle redshift distortions agrees with simulations, and that 
it represents an improvement over previous works, at modest 
cost in compexity . Poss ible generalizations along the lines of 
iMatsubara et all (|2004l ) are left for future research. 
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